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Rapidity and transverse momentum dependence of inclusive J/i/a 
production in pp collisions at TeV 



The ALICE experiment at the LHC has studied inclusive J/y/ production at central and forward 
rapidities in pp collisions at ^Js = 7 TeV. In this Letter, we report on the first results obtained detecting 
the J/y/ through the dilepton decay into e+e^ and jJ-^jJ.^ pairs in the rapidity ranges \y\ < 0.9 and 
2.5 < y < 4, respectively, and with acceptance down to zero pj. In the dielectron channel the analysis 
was carried out on a data sample corresponding to an integrated luminosity Lin, = 5.6 nb^^ and the 
number of signal events is Nj/^ = 352 ± 32 (stat.) ± 28 (syst.); the corresponding figures in the 
dimuon channel are Lin, = 15.6 nb^' and Nj/y, = 1924 ± 77 (stat.) ± 144 (syst.). The measured 
production cross sections are CJj/^^dyl < 0.9) = 12.4 ±1.1 (stat.) ± 1.8 (syst.) (syst.pol.) /lb 
and C7j/,^(2.5 < y < 4) = 6.31 ± 0.25 (stat.) ± 0.76 (syst.) t^lH (syst.pol.) jUb. The differential cross 
sections, in transverse momentum and rapidity, of the J/\// were also measured. 
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1 Introduction 

The hadroproduction of heavy quarkonium states is governed by both perturbative and non-perturbative 
aspects of Quantum Chromodynamics (QCD) and was extensively studied at the Tevatron [T]|2]|3]lll and 
RHIC [5J hadron colliders. Various theoretical approaches, recently reviewed in ||6]|3, were proposed 
to describe the data. They mainly differ in the details of the non-perturbative evolution of the heavy 
quark pair towards a bound state. The models are not able to consistently reproduce the production cross 
section, the transverse momentum (pj) distributions and the polarization. Recently, theoretical studies 
focused on the calculation of NLO and NNLO contributions, finding that their impact on the results is 
quantitatively important ||8] |9] [10] [11] |T211. Measurements in the new energy domain of the LHC are 
crucial for a deeper understanding of the physics involved in hadroproduction processes. Furthermore, 
the range of Bjorken-x values accessible at LHC energies is unique. Low-pj charmonium measurements, 
in particular at forward rapidity, are sensitive to an unexplored region (x <10^^ at = >n^/^) of the 
gluon distribution function of the proton. 

Heavy quarkonia are measured in the ALICE experiment |[T3l through their e+e^ and pL^pi^ decays. In 
this Letter we present the results on inclusive J/i//^ production in pp collisions at = 7 TeV, measured 
in the rapidity regions \y\ < 0.9 for the dielectron channel and 2.5 < y < 4.0 for the dimuon decay 
channeU. First, a description of the ALICE experimental apparatus is given, mainly focusing on the muon 
detection. Track reconstruction in the central rapidity region was discussed previously |[T4l . Details are 
provided concerning the data analysis, the reconstruction algorithm, the event selection criteria and the 
techniques used for the extraction of the signal. After describing the determination of the acceptance 
and efficiency corrections, and the methods used for the evaluation of the luminosity, the values of the 
integrated, ^-differential and /^x-differential J/y/ cross sections are presented and compared with the 
results obtained by the other LHC experiments |[T5][T6][r7l . 

2 Experimental apparatus and data taking conditions 

The ALICE experiment |[T3l consists of two main parts: a central barrel and a muon spectrometer. The 
central barrel detectors (|t]| <0.9) are embedded in a large solenoidal magnet, providing a magnetic 
field of 0.5 T. Various detector systems track particles down to pj of about 100 MeV/c, and can provide 
particle identification over a wide momentum range. The muon spectrometer covers the pseudorapidity 
range — 4< rj < —2.5 and detects muons having p >4 GeV/c. Finally, various sets of forward detectors 
further extend the charged particle pseudorapidity coverage up to T] = 5.1 and can be used for triggering 
purposes. 

The barrel detectors used in this analysis are the Inner Tracking System (ITS) and the Time Projection 
Chamber (TPC). The ITS [.13. 181 is a cylindrically-shaped silicon tracker that surrounds the central 
beam pipe. It consists of six layers, with radii between 3.9 cm and 43.0 cm, covering the pseudo-rapidity 
range |tj| < 0.9. The two innermost layers are equipped with Silicon Pixel Detectors (SPD), the two 
intermediate layers contain Silicon Drift Detectors (SDD), and Silicon Strip Detectors (SSD) are used on 
the two outemiost layers. The main task of the ITS is to provide precise track and vertex reconstruction 
close to the interaction point, to improve the overall momentum resolution and to extend tracking down 
to very low pj. The SPD can also deliver a signal for the first level trigger (LO), which is based on a hit 
pattern recognition system at the level of individual readout chips. 

The TPC |jT9ll is a large cylindrical drift detector with a central high voltage membrane maintained at 
— 100 kV and two readout planes at the end-caps. The active volume extends over the ranges 85 < r < 
247 cm and —250 < z < 250 cm in the radial and longitudinal (beam) directions, respectively. Besides 

'The muon spectrometer covers, in the ALICE official reference frame, a negative T] range and, consequently, a negative 
y range. However, since in pp the physics is symmetric with respect to y=0, we have dropped the negative sign when quoting 
rapidity values. 
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being the main tracking detector in the central barrel, the TPC also provides charged hadron identification 
with very good purity up to a total momentum of about 3 GeV/c ll20l and electrons up to about 10 
GeV/c, via specific energy loss (dE/dx) measurement, with a resolution of a=5.5% for minimum ionizing 
particles. 

Other central barrel detectors with full azimuthal coverage, in particular the Time-Of-Flight (TOP) ETTl . 
the Transition Radiation Detector (TRD) [22] and the Electromagnetic Calorimeter (EMCAL) [23], sue 
not used in this analysis, but are expected to significantly improve the electron identification and trigger- 
ing capabilities of the experiment in the future. 

The muon spectrometer consists of a front absorber followed by a 3 T-m dipole magnet, coupled to 
tracking and triggering devices. Muons emitted in the forward rapidity region are filtered by means of a 
10 interaction length (Ai) thick front absorber made of carbon, concrete and steel, and placed between 0.9 
and 5.0 m from the nominal position of the interaction point (IP). Muon tracking is performed by means 
of 5 tracking stations, positioned between 5.2 and 14.4 m from the IP, each one based on two planes of 
Cathode Pad Chambers. The total number of electronic channels is close to 1.1 • 10^, and the intrinsic 
spatial resolution for these detectors is of the order of 70 jxra in the bending direction. Stations 1 and 2 
(4 and 5) are located upstream (downstream) of the dipole magnet, while station 3 is embedded inside its 
gap. A muon triggering system is placed downstream of a 1.2 m thick iron wall (7.2 Ai), which absorbs 
secondary hadrons escaping the front absorber and low-momentum muons (having p < 1.5 GeV/c at the 
exit of the front absorber). It consists of two stations positioned at 16.1 and 17.1 m from the IP, equipped 
with two planes of Resistive Plate Chambers (RPC) each. The spatial resolution achieved is better than 
1 cm, while the time resolution is of the order of 2 ns. Throughout its full length, a conical absorber 
{6 < 2°) made of tungsten, lead and steel protects the muon spectrometer against secondary particles 
produced by the interaction of large- tj primaries in the beam pipe. 

Finally, the VZERO detector consists of two scintillator arrays covering the range 2.8 < tj < 5.1 and 
— 3.7 < T] < — 1.7, and positioned, respectively, at z = 340 and z = —90 cm from the IP. It provides 
timing information to the LO trigger with a resolution better than 1 ns. This feature proves to be useful 
in the offline rejection of beam-halo and beam-gas events. 

The results presented in this Letter were obtained by analyzing data collected in the first year of operation 
of the LHC, corresponding to pp collisions at ^/s=l TeV. During this period, the LHC reached its goal of 
delivering more than 10^^ cm^^s^^ instantaneous luminosity. In ALICE, the instantaneous luminosity 
was kept to 0.6 — 1.2 • 10^^ cm^^ s^^ in order to have a collision pile-up rate in the same bunch crossing 
below 5%. In order to increase the statistics for low cross-section processes, ALICE ran, during short 
periods, at luminosities about 10 times higher, therefore with a much larger pile-up rate in the same bunch 
crossing. For the analysis presented in this Letter, the collected data were divided into three sub-periods. 
Each one corresponds to similar average luminosity and pile-up rates, and is characterized by reasonably 
stable tracking and trigger detector configuration and performance. 

The event sample used in this analysis corresponds to minimum bias events (MB trigger) and, for the 
muon analysis, to events where the detection of at least one muon in the angular acceptance of the 
muon spectrometer (jX-MB trigger) is additionally required. The MB trigger is defined as the logical 
OR between the requirement of at least one fired readout chip in the SPD, and a signal in at least one 
of the two VZERO detectors [24|. It also requires a coincidence with signals from two beam pick-up 
counters, one on each side of the interaction region, indicating the passage of proton bunches. The jj.- 
MB trigger allows the selection of events where at least one particle was detected in the trigger chambers 
of the muon spectrometer. The trigger logic is based on the requirement of having at least 3 hits (out 
of 4) in the trigger stations, both in the bending and non-bending directions [25 |. In this way one can 
define a "trigger track", and compute its deviation with respect to a track with infinite momentum. By 
requiring such a deviation to be smaller than a certain value one can select muon candidate tracks having 



Rapidity and transverse momentum dependence of inclusive J/y production in pp . 



5 



a transverse momentum larger than a pre-defined value. Such a Pj. cut can be used to reject soft muons, 
dominated by 7i and K decays, and is able to limit the muon trigger rate when the machine luminosity 
is high. The instantaneous luminosity at ALICE allowed for the choice of the lowest threshold 
(0.5 GeV/c), leading to a /x-MB trigger rate between 30 and 500 Hz. With this pj^, the effect of the 
trigger response function on the J/i/a detection efficiency is negligible. Note that the effect of such a cut 
is not sharp and that the selection efficiency reaches the plateau value only at pj GeV/c. Finally, in 
order to limit the systematic uncertainties related to non- uniformities in the detector response, data were 
collected by periodically varying the polarities of the solenoidal and dipole magnets. 

3 Data analysis 

For the dielectron analysis, 3.5-10^ minimum bias events (A^mb) are analyzed. An event with a recon- 
structed vertex position Zy is accepted if |zv| < 10 cm. The tracks are required to have a minimum pj of 
1.0 GeV/c, a minimum number of 70 TPC clusters per track (out of a maximum of 159), a per space 
point of the momentum fit lower than 4, and to point back to the interaction vertex within 1 cm in the 
transverse plane. A hit in at least one of the two innermost layers of the ITS is required in order to reduce 
the contribution of electrons from 7 conversions. For full-length tracks, the geometrical coverage of the 
central barrel detectors is |t]| < 0.9. 




20 1 < — < < < — 1 

0.3 0.4 1 2 3 4 5 6 7 8 91 20 

p (GeV/c) 



Fig. 1: Specific energy loss in the TPC as a function of momentum with superimposed Bethe-Bloch lines for 
various particle species. The dashed lines show the pion and proton exclusion bands. The dotted line corresponds 
to the +3a cut for electrons (see text). 

The particle identification performance of the TPC is essential for the J/\j/ measurement. In Fig. [1] the 
specific energy loss in the TPC is shown as a function of momentum in the region of interest for the 
present measurement. A ±3 a inclusion cut for electrons and a ±3.5 a (±3 a) exclusion cuts for 
pions (protons) were employed. As seen in Fig. [T] with our current identification strategy, the electron 
identification is performed with an efficiency better than 50% for momenta below 7-8 GeV/c. 

Electron candidates compatible, together with a positron candidate, with being products of 7 conversions 
were removed, in order to reduce the combinatorial background. It was verified, using a Monte Carlo 
simulation, that this procedure does not affect the J/i/A signal. 

The invariant mass distribution for the opposite-sign (OS) electron pairs is shown in Fig. |2] In the 
same figure we also show the background contribution, obtained as the sum of the like-sign (LS) pairs, 
N^^ +N , scaled to match the integral of the OS distribution in the mass interval 3.2—5.0 GeV/c^. 
The scale factor, 1.23, originates from the presence of correlated background (mostly from semilep- 
tonic charm decays) in the OS distribution, but is also influenced by misidentified electrons and by 
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Fig. 2: Top panel: invariant mass distributions for opposite-sign (OS) and like-sign (LS) electron pairs {\y\ < 0.9, 
all pj), as well as for pairs obtained with one track randomly rotated (TrkRot, see text). Bottom panel: the 
difference of the OS and LS distributions with the fit to the Monte Carlo (MC) signal superimposed. 

electrons from conversions. In the top panel of Fig. |2] we also show the background estimated using 
a track rotation method (TrkRot) H, used later in the estimate of the systematic uncertainties related to 
signal extraction. The signal, obtained by subtracting the scaled LS distribution from the OS, is shown 
in the bottom panel of Fig. |2]in comparison with the signal from Monte Carlo (MC) simulations (de- 
scribed below). A good agreement between data and MC is observed, both for the bulk of the signal 
and for the bremsstrahlung tail. Integration of the signal in the mass range 2.92—3.16 GeV/c^ yields 
Ny^ - 352 lb 32 (stat.) it 28 (syst.) counts (the systematic uncertainty on this quantity is described 
below); the signal to background ratio is S/B = 1.2 it 0.1 and the significance is S/\/S + B = 13.9 ± 0.6. 
The taggingand corresponding rejection of / conversions is found to improve S/B by ~30%. The MC 
simulations □ show that (66.8 it 1.9)% of the signal is within the integration range. The main contribution 
to the uncertainty on this quantity was obtained by analyzing MC samples where the detector material 
budget was varied by ±6% ll29l with respect to the nominal value, and by varying the track-related cuts 
(Pj and required number of TPC clusters) around their nominal values. A smaller contribution (1 %, in 
terms of the relative error) is attributed to the small discrepancies between the invariant mass distribution 
as provided by QED at the next to leading order [30.1 and by the event generator (EvtGen -i- PHOTOS); 
the latter contribution remains even after taking into account the detector resolution. A fit to the invariant 
mass distribution after background subtraction with a Crystal Ball function lISTl gives a mass resolution 
of 28.3±1.8MeV/cl 

For the dimuon channel, the total data sample available for physics analysis amounts to 1.9-10^ MB 
events, of which 1.0-10^ satisfy the ^u-MB condition. 

^The method consists in rotating, around the z axis, one of the tracks of the OS pair by a random azimuthal angle. More 
pairs can be obtained by applying the method several times to the same pair. 

■'in the MC simulation the decay of the particles was handled by the EvtGen package 1261 and the final state radiation 
was described using PHOTOS I27||28| . 



Rapidity and transverse momentum dependence of inclusive J/y production in pp . 



7 



An accurate alignment of the tracking chambers of the muon spectrometer is an essential pre-requisite to 
identify resonances in the pL^pi^ invariant mass spectrum. This was carried out using a modified version 
of the MILLEPEDE package ll32l[33l . starting from a sample of 3-10^ tracks, taken with no magnetic 
field in the dipole and in the solenoid. The resulting alignment precision is ~750 /im in the bending and 
non-bending directions. 

Track reconstruction is based on a Kalman filter algorithm ll33l[34l . The procedure starts from the most 
downstream tracking stations (4 and 5), which are less subject to the background due to soft particles 
that escape the front absorber. Straight line segments are formed by joining clusters on the two planes of 
each station and a first estimate of the track parameters (position, slope and inverse bending momentum) 
and corresponding errors is made. The momentum is first estimated assuming that the track originates 
from the vertex and is bent by a constant magnetic field in the dipole. In a second step, track candidates 
on station 4 are extrapolated to station 5 (or vice versa) and paired with at least one cluster on the basis 
of a cut. If several clusters are found, the track is duplicated to consider all the possible combinations. 
After this association the track parameters and errors are recalculated using the Kalman filter. 

The same procedure is repeated iteratively for the upstream stations, rejecting, at each step, the candidates 
for which no cluster is found or those whose parameters indicate that they will exit the geometrical 
acceptance of the spectrometer in the next steps. At the end of the procedure, additional algorithms are 
applied to improve the track quality by adding/removing clusters based on a cut, and removing fake 
tracks sharing clusters with others. Finally, the remaining tracks are extrapolated to the primary vertex 
position as given by the SPD |24], and their parameters are recomputed taking into account the energy 
loss and multiple Coulomb scattering in the absorber. With the alignment precision obtained for the 
analyzed data sample the relative momentum resolution of the reconstructed tracks ranges between 2% 
at 10 GeV/c and 10% at 100 GeV/c. 

After reconstruction, 4.1-10^ events having at least two muon candidates are found, out of which only 
6% have three or more muons. Various selection cuts are then applied to this data sample. First, events 
are required to have at least one interaction vertex reconstructed by the SPD. This cut rejects 0.5% of the 
statistics. Then, it is required that at least one of the two muon candidates matches the corresponding 
hits in the trigger chambers. In this way hadrons produced in the absorber, which are stopped by the iron 
wall positioned upstream of the trigger chambers, are efficiently rejected. This cut rejects ~24% of the 
muon pairs, and its effect is important only for m^^ <1 GeV/c^. In fact, since in 99% of the cases at least 
one of the two J/i/a decay muons has a transverse momentum larger than the trigger /jt threshold, the 
signal loss induced by this cut is negligible. Requiring both candidate muon tracks to be matched with 
the corresponding "trigger tracks" would increase the purity of the muon sample, but it was checked 
that this cut would lead to a loss of ~ 20% of the J/i/a events without decisively increasing the signal 
to background ratio at the J/i/A invariant mass. Furthermore, the cut /?abs > 17.5 cm, where /?abs is the 
radial coordinate of the track at the end of the front absorber, was applied. In this way, muons emitted 
at small angles, that have crossed a significant fraction of the thick beam shield, can be rejected. Finally, 
to remove events very close to the edge of the spectrometer acceptance, the cut 2.5 < 3^ < 4 on the pair 
rapidity was applied. These quality cuts reject 10.3% of the muon pairs. 

After selection, the dimuon sample consists of 1.75-10^ OS muon pairs. In Fig. |3]we present the OS 
invariant mass spectrum for the mass region 1.5 < m^^ < 5 GeV/c^, corresponding to the sub-period 
having the largest statistics. A peak corresponding to the J/y M^M decay is clearly visible in the 
spectrum, on top of a large continuum. A weaker signal, corresponding to the '//^(25) decay, is also 
visible, in spite of the poor signal to background ratio. 

The number of signal events A'^j was extracted by fitting the mass range 1 .5 < m^^ < 5 GeV/cl The J/i/a 
and V'^(25) line shapes are described with Crystal Ball functions [3iJ, while the underlying continuum 
was parameterized using the sum of two exponentials. The functions representing the resonances were 
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Fig. 3: Invariant mass distribution for opposite-sign muon pairs (2.5 < y < 4, all pj), in the mass region 1.5 < 
wi/j/j < 5 GeV/c^ with the result of the fit. The plot refers to the sub-period with the largest statistics (A'j/^ = 
957 ±56, corresponding to Lint = 7.9 nb^'). The fitted J/y/ and i/(25) contributions, as well as the background, 
are also shown. 

obtained by fitting the expected mass distribution of a pure MC signal sample. Such a sample was 
obtained by generating, for each sub-period, 3/\j/ and ^^(25) events with realistic differential distributions 
(see below for details). In order to account for small uncertainties in the MC description of the set-up, 
the position of the I/i/a mass pole, as well as the width of the Crystal Ball function, were kept as free 
parameters in the invariant mass fit. Due to the small statistics, the '//^(25) parameters were tied to those 
of the J/y, imposing the mass difference between the two states to be equal to the one given by the 
Particle Data Group (PDG) |[35l . and the ratio of the resonance widths to be equal to the one obtained in 
the MC. 

This choice of parameters leads to a satisfactory fit of the invariant mass spectrum (X'^ /ndf =1.14), as 
shown in Fig.[3j The Crystal Ball function describing the J/xj/ is peaked at /mj/,^ = 3.1 18 it 0.005 GeV/c^. 
Such a value is larger than the one quoted by the PDG group by only 0.6%, showing that the accuracy of 
the magnetic field mapping and of the energy loss correction is reasonably under control. The measured 
width of the Crystal Ball function is aj/^, = 94 it 8 MeV/c^, in agreement within less than 2% with the 
MC, and its FWHM is 221 MeV/c^. 

The same fitting procedure, applied to the other sub-periods, gives consistent results in terms of both mj 
(within 0.2%) and Oy^, (within 4%). The signal to background ratio, in the mass range 2.9 < m^^ < 3.3 
GeV/c^, varies between 2.3 and 2.9 in the various sub-periods. The total number of J/y signal events, 
obtained by integrating the Crystall Ball function over the full mass range, is Nj^^ = 1924 it 77(stat.) it 
144(syst.). The determination of the systematic uncertainty on A'j/i^ is described later in Section|5] 

4 Acceptance and efficiency corrections, luminosity normalization 

In order to extract the J/\j/ yield, the number of signal events must be corrected, with a MC procedure, 
for the acceptance of the apparatus and for reconstruction and triggering efficiencies. This procedure 
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is based on the generation of a large sample of signal events, with a pj distribution extrapolated from 
CDF measurements [ 1 1 and a y distribution parameterized from Color Evaporation Model (CEM) calcu- 
lations Il36l . To avoid the loss of events due to smearing effects at the edge of the angular acceptance, the 
generation was performed over y ranges wider than those covered by the two detector systems. It was 
also assumed that J/y production is unpolarized. The acceptance factors are obtained with respect to the 
J/Y rapidity ranges |y| < 0.9 and 2.5 < y < 4.0 for the central barrel and muon detectors, respectively. 

For the central barrel detectors, the acceptance times efficiency value (A x e) is 8.5% and is the product 
of four contributions: i) a kinematic factor, namely the requirement of having both e+ and e^ within 
the acceptance (|t]'^ | < 0.9), satisfying a transverse momentum cut ''^ > 1 GeV/c. This factor 
amounts to 31.3%; ii) the reconstruction efficiency for the e+e^ pair, which is 50.3%; iii) the identifica- 
tion efficiency, which is 81.0%; iv) the fraction of the signal within the mass range 2.92-3.16 GeV/c^, 
which is 66.8%. 

For the muon spectrometer, the tracking efficiency is calculated from MC simulations, including a re- 
alistic map of dead channels and the residual misalignment of the detection elements. This efficiency, 
obtained from a sample of generated tracks which match the hit pattern on the various chambers required 
by the reconstruction algorithm, is (97.1 ±0.8)%. 

The efficiencies of the muon trigger chambers are obtained from the analysis of the "trigger tracks" 
collected in the measured data sample. The "trigger tracks", as explained in Section 2, are defined by 
the presence of a hit in at least 3 (out of 4) trigger planes. The efficiency for a chamber belonging to 
a certain trigger plane is calculated starting from a sample of "trigger tracks" where the corresponding 
chambers on the other 3 planes have recorded a hit, and then looking for the presence or absence of 
a hit in the chamber under study. Such a requirement on the sample does not introduce a bias in the 
efficiency calculation because the response of the detector planes are independent. Typical efficiency 
values are around 96%, with 90% of the detector surface having an efficiency larger than 91%. The 
obtained efficiencies are then plugged in the simulations in order to provide a realistic description of 
the detector. The time variation of the tracking and trigger detector efficiencies was accounted for in 
the simulation. Internally to each sub-period, the response of the tracking chamber channels is further 
weighted run by run. 

For each sub-period /, the ratio between the total number of reconstructed events, satisfying the anal- 
ysis cuts, and the generated events in the range 2.5 < y <4 gives the product A x Si for the J/y/. The 
differences between the Ax Si for the various sub-periods do not exceed 8%, and their average value is 
{Ax e) = 32.9%. It is worth noting that Axe exhibits, for both the central barrel detectors and the muon 
spectrometer, a rather small variation as a function of the J/i// pj, down to zero pj. 

To get the production cross section value, the ratio N^J^^ = Ny^/{A x e) must be normalized to the 
integrated luminosity, or to the measured cross section for a chosen reference process. For this analysis, 
the adopted reference is the occurrence of the MB condition itself. One has simply 

BR{J/Y ^ in-) Nub ^ ^ 

where BR(J/i/a — £+£^)=(5.94 it 0.06)% |[35l . A^mb is the number of minimum bias collisions and Gmb 
is the measured cross section for such events. A^mb was corrected, run by run, for the probability of 
having multiple interactions in a single bunch crossing. 

In the muon channel, the J/i/a signal was collected using the jX-MB trigger condition. Therefore, Eq. [T] 
has to include a multiplicative factor R that links the occurrence of a reference process in the jX-MB and 
MB event samples. We have chosen as a reference process the yield A'^^ of single muons, detected in the 
region -4 < T] < -2.5, and with px > 1 GeV/c. The R factor is then defined as the ratio R = N^^/N^^^^^ 
of the single muon yields for the two event samples. The numerical values of the R factor strongly depend 
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on the relative bandwidth assigned by the data acquisition to the two trigger samples in each sub-period, 
and vary between 0. 10 and 0.42. However, the choice of the pj cut has no significant influence on these 
values (<1% for cut values between and 3 GeV/c). This is due to the fact that both the jX-MB and 
MB muon samples are subject to the same set of cuts, including the requirement of matching of the 
reconstructed track with the corresponding trigger track. 

The Gmb value is 62.3 mb, and is affected by a 4% systematic uncertainty. It was obtained relative to the 
cross section avoAND f3T|, measured in a van der Meer scan [,38 j . of the coincidence VOAND between 
signals in the two VZERO detectors. The relative factor cTvoand /omb was obtained as the fraction of 
MB events where the LO trigger input corresponding to the VOAND condition has fired. Its value is 0.87, 
and is stable within 0.5% over the analyzed data sample. The integrated luminosity Li^ = A^mb/cmb is 
5.6 nb^' for the dielectron sample. For the dimuon sample Ljnt = (A^mb/o'mb)/^ = 15.6 nb~^ 

5 Systematic uncertainties 

The systematic uncertainty on the inclusive J/y cross section measurement was obtained considering the 
following sources: 

- The uncertainty on the signal extraction procedure, for the electron channel, was estimated using 
the track rotation method as an alternative background calculation procedure, see Fig. |2j and by fit- 
ting the invariant mass distributions with a convolution of a polynomial and a Crystal Ball function 
(with parameters constrained via Monte Carlo). We have also varied the invariant mass ranges for 
the signal extraction and for background normalization. The value obtained is 8%. Including the 
contribution from the uncertainty on the material budget and the accuracy of the event generator 
in describing the radiative decay of the J/i/a (J/y — )-e+e^7) leads to a value of 8.5%. For the muon 
channel, various tests were performed. In particular, we tried to release in the fit the values of the 
parameters governing the asymmetric left tail of the Crystal Ball line shape, which were fixed to 
their MC values in the default fitting procedure. Alternative functions for the description of the 
signal and background shapes were also used. In particular, for the J/y a variable-width Gaussian 
function, adopted in the past by the NA50 and NA60 Collaborations [39], was used. For the back- 
ground, a different shape was tested, based on a Gaussian having a width continuously increasing 
with the mass. The estimated overall systematic uncertainty on the signal extraction is 7.5%. 

- The acceptance calculation depends on the y and pj input distributions. For the electrons, the 
uncertainty is mainly determined by the choice of the pj spectrum. By varying the {pj) of the input 
distribution within a factor 2, a 1.5% variation in the acceptance was obtained. Such a small value 
is indeed a consequence of the weak pj dependence of the acceptance for the bulk of the spectrum. 
For the muons, both y and pj were varied, using as alternative distributions those expected for pp 
collisions at ^/s=4 and 10 TeV BOl . As a further test, the measured differential distributions 
obtained from the analysis described later in Section 6, were used as an input in the calculation 
of the acceptance. In this way, a 5% systematic uncertainty on this quantity was determined. The 
larger systematic uncertainty for the muon channel is due to the larger influence, for this channel, 
of the choice of the shape of the rapidity distribution. 

- The uncertainty on the muon trigger efficiency calculation was estimated comparing NjJ^^ for the 
sample where only one of the two decay muons is required to match the trigger condition, with 
the same quantity for the sample where both muons satisfy that condition. The 4% discrepancy 
between the two quantities is taken as the systematic uncertainty on the evaluation of the trigger 
efficiency. 

- The uncertainty on the reconstruction efficiency, for the central barrel analysis, is due to the track 
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quality (4%) and particle identification (10%) cuts and originates from residual mismatches be- 
tween data and MC simulations. 

For the muon analysis, the systematic uncertainty can be estimated by comparing determinations 
of the tracking efficiency based on real data and on a MC approach. In the first case, the tracking 
efficiency can be evaluated starting from the determination of the efficiency per chamber, com- 
puted using the redundancy of the tracking information in each station. The values thus obtained 
are in the range from 91.8 to 99.8%. The tracking efficiency evaluated starting from these chamber 
efficiencies is (98.8 ±0.8)%. The very same procedure, in a MC approach, gives a (99.8^q'3)% 
tracking efficiency. These two quantities differ by 1%, which is taken as an estimate of the system- 
atic uncertainty. However, this method is not able to detect losses of efficiency due to the presence 
of correlated dead-areas in the same region of two chambers belonging to the same station. Such 
correlated dead-areas were singled out by studying, on data, the cluster maps of each station, and 
the corresponding loss of efficiency was estimated to be (2.8 ±0.4)%. Taking into account this 
effect, the resulting tracking efficiency is in good agreement with the value previously quoted (see 
Section|4ll from realistic MC simulations, (97.1 ±0.8)%. Nevertheless, a 1% additional system- 
atic error (30% of the efficiency loss discussed above) was assumed, to take into account possible 
small-area correlations that could be missed in the present approach. Combining this error with 
the one previously mentioned, the overall systematic uncertainty on the muon tracking efficiency 
is 1.5%, which gives 3% for muon pair detection. 

- The error on the luminosity measurement is dominated by the 4% systematic uncertainty on the 
determination of avoAND> which is due to the uncertainties on the beam intensities PTI and on the 
analysis procedure related to the van der Meer scan of the VOAND signal. Other effects, such as 
the oscillation in the ratio between the MB and VOAND counts, contribute to less than 1%. The 
cross section o^^mb, relative to the occurence of the /i-MB trigger, was also measured in a van 
der Meer scan Il37ll38l and was used as an alternative reference for the luminosity determination 
in the muon analysis. Using a^-MB as a reference cross section, a 4% difference has been found 
with respect to the integrated luminosity based on ajviB- For safety, this 4% has been quadratically 
added to the luminosity systematic error in the muon analysis. In addition, for the muons, the 
calculation of the integrated luminosity, as described above, is also connected with the estimate 
of the R factor. This quantity was evaluated in an alternative way, using the information from the 
trigger scalers and taking into account the dead-time of the triggers. By comparing the two results, 
a 3% systematic uncertainty on the R factor was estimated. 

- The branching ratio of the J/i// decay to lepton pairs is known with a 1% accuracy. 

- The acceptance values significantly depend on the degree of polarization assumed in the J/\j/ distri- 
butions. They were calculated in the two cases of fully transverse (A = 1) or longitudinal (A = — 1) 



polarizatioqj, in the CoUins-Soper (CS) and helicity (HE) reference frames. 

The systematic uncertainties are summarized in Table [T] The systematic uncertainty on the inclusive J/i/a 
cross section is obtained by quadratically combining the errors from the sources described above, except 
polarization, and is 12.1% for the dimuon channel and 14.5% for the dielectron one. The systematic 
uncertainty due to the unknown J/y polarization will be quoted separately. 

6 Integrated and differential J/xj/ cross sections 



The inclusive J/i//^ production cross sections in pp collisions at ^/s=^ TeV are: 

(^}/xifi\y\ < 0.9) = 12.4 ± 1.1 (Stat.) ± 1.8 (syst.) ± 1.8 (Ahe = 1) - 2.7 (Ahe = -1) Mb and 

c7j/^(2.5 < 3; < 4) = 6.31 ± 0.25 (stat.) ± 0.76 (syst.) + 0.95 (Acs = 1) - 1-96 (Acs = - 1) jUb- 




'^The polar angle distribution of the decay leptons is given by dN/d cos 9 — l+X cos^ 6 
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Table 1: Systematic uncertainties (in percent) on the quantities associated to the integrated J/y/ cross section 
measurement. 



Channel 


e+e 




Signal extraction 


8.5 


7.5 


Acceptance input 


1.5 


5 


Trigger efficiency 





4 


Reconstruction efficiency 


11 


3 


R factor 




3 


Luminosity 


4 


5.5 


B.R. 


1 


Polarization 


A = -1 


A = l 


A = -l 


A = 1 


CS 


+19 


-13 


+31 


-15 


HE 


+21 


-15 


+22 


-10 




mee (GeV/c^) 



Fig. 4: Invariant mass spectra for OS electron pairs ( \y\ < 0.9), in bins of pj. The background calculated using LS 
and TrkRot approaches are also shown. 
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Fig. 5: Invariant mass spectra for OS muon pairs (2 .5 < y < 4), in bins of pj. The results of the fits are also shown. 

The systematic uncertainties related to the unknown polarization are quoted for the reference frame 
where they are larger. 

In the dielectron channel, the dOj/^/dp-Y differential cross section was measured in five pj bins, between 
and 7 GeV/c. In each bin, the signal was extracted with the same approach used for the integrated 
invariant mass spectrum. In Fig. |4]the OS invariant mass spectra are shown, together with the LS and the 
TrkRot backgrounds. The corrections for acceptance and reconstruction efficiency and the systematic 
errors are given in Table |2] Some of the contributions to the systematic uncertainty do not depend on pj, 
thus affecting only the overall normalization, and they are separately quoted in Table|2] The contributions 
which depend on pj, even when they are correlated bin by bin, were included among the non-correlated 
systematic errors. 

For the analysis in the dimuon channel, a differential study of J/y production was performed in the two 
kinematic variables y and pj separately. In particular, daj/^/d/^x was studied in seven bins between 
and 8 GeV/c, and dOj/^/dy in five bins between 2.5 and 4. The event sample used for the determination 
of the differential cross sections is slightly smaller (by about 15%, corresponding to Lint=13.3 nb^^) than 
the one analyzed for the integrated cross section. This is due to the fact that the statistics in one of the 
three sub-periods of the data taking is too small to allow a satisfactory fit of the differential invariant 
mass spectra. 

The J/y signal was extracted, for each y or pj bin, with the same fitting technique used for the integrated 
invariant mass spectra. Since the V'^(25') yield is rather small and cannot be safely constrained by the data 
themselves, its contribution was fixed in such a way as to have the same va(25')/(J/i//^) ratio extracted 
from the integrated spectrum. Anyway, the results of the fit, for what concerns Ny^, are quite insensitive 
to the precise level of the '/'^(25) contribution. It has been verified, for example, that fixing for each pj 
bin the ratios V'^(25') /{3/y) to the values measured (in the range pt >2 GeV/c) by CDF ll42]| . A'j/i^ varies 
by less than 1%. In Fig.[5]the OS invariant mass spectra are shown, together with the result of the fits. 

The acceptance times efficiency was calculated differentially in y and pj and the values are reported in 
Table 121 It can be noted that as a function of pj, the Axe coverage of the muon spectrometer for J/i/a 
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Fig. 6: d^ffj/i^/dpxdy for the midrapidity range and for the forward rapidity data, compared with results from the 
other LHC experiments ifT5] [T6l [TTl . obtained in similar rapidity ranges. The error bars represent the quadratic 
sum of the statistical and systematic errors, while the systematic uncertainties on luminosity are shown as boxes. 
The symbols are plotted at the center of each bin. 



production extends down to zero pj, and that the values vary by less than a factor 1.6 in the analyzed pj 
range. Axe has a stronger y dependence, but its values are larger than 10% everywhere. 

The differential cross sections are then calculated with the same approach used for the integrated cross 
section, normalizing Nj°'^{y) and NjJ^^{pt) to the integrated luminosity. The differential cross sections 
are affected by the same systematic error sources discussed in the previous section. All except the one 
related to the signal extraction can be considered as common, or strongly correlated. Table |2] gives a 
summary of the results, including the various sources of systematic uncertainties (correlated, uncorre- 
cted and polarization-related). 

The results are presented in Fig. |6]and Fig. |7]for the /^T-differential cross section (f-ay^i,/dpjdy and 
dCj/vz/dj (pt > 0), respectively. For the rapidity distribution, the values obtained in the forward region 
were also reflected with respect to j = 0. In both figures, the symbols are plotted at the center of each 
bin. The statistical and systematic errors were added in quadrature, apart from the 4 (5.5)% systematic 
uncertainties on luminosity for the dielectron (dimuon) channels, shown as boxes. The differential cross 
sections shown in Fig. [6] and |7] assume unpolarized J/y production. Systematic uncertainties due to 
the unknown J/xj/ polarization are not shown. Our results are compared with those by the CMS 1 15]]. 
LHCb 1161 and ATLAS 1 171 Collaborations. Also for these data the uncertainty due to luminosity, which 
is 1 1% for CMS, 10% for LHCb and 3.4% for ATLAS, is shown separately (boxes), while the error bars 
contain the statistical and the other sources of systematic errors added in quadrature. Our measurement at 
central rapidity reaches pj = and is therefore complementary to the data of CMS, available at |y| < 1.2 
for Pj > 6.5 GeV/c, and ATLAS, which covers the region \y\ < 0.75, pj >7 GeV/c. In order to compare 
our d^Gy^/dpTdy in the forward rapidity range with that of LHCb, we added the LHCb data for prompt 
and non-prompt production and integrated in the range 2.5 < j < 4 to match our measurement. The 
agreement between the two data sets is good. 

In Fig. |7] our results are compared with the corresponding values from the CMS and LHCb experiments. 
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Fig. 7: dGj/^/dy, compared with results from the other LHC experiments llT5][T6l[T7l . The error bars represent the 
quadratic sum of the statistical and systematic errors, while the systematic uncertainties on luminosity are shown 
as boxes. The symbols are plotted at the center of each bin. 

for the rapidity bins where the pj coverage extends down to zero (ATLAS has no coverage down to 
Pj=0 in any rapidity range). For CMS, the value for 1.6 < \y\ < 2.4 was obtained by integrating the 
pubhshed d^cyj/^/dpxdj data ifTSl . while for LHCb the pubhshed doy^/dy for prompt and non-prompt 
production fl^ were added. Our data, together with that of the other LHC experiments, constitute a 
comprehensive measurement of inclusive J/y production cross section as a function of rapidity. At the 
LHC, the inclusive J/i/A production cross section at central rapidity is almost twice larger than at Tevatron 
(V?=1.96 TeV) m and about ten times larger than at RHIC (^^=0.2 TeV) 10. The width (FWHM) of 
the rapidity distribution derived from our data is about twice larger than at RHIC Q. 

We stress that the results described in this Letter refer to inclusive J/y production. Therefore the mea- 
sured yield is a superposition of a direct component and of J/y coming from the decay of higher-mass 
charmonium states, in particular the Xci, Xc2 and Yi^^) states. These contributions were measured in 
lower-energy experiments and were found to be ~25% (Xci + Xci) and ~8% (yC^'^)) of the total mea- 
sured J/y yield B3l |44|| . The XcO contribution is negligible since its B.R. into 3/y is of the order of 
1%. In addition to this "prompt" production, decays of beauty hadrons are also known to give a sizeable 
contribution (of the order of 10-15% in the pj range accessed by ALICE [16]) to the observed J/y yield. 
With future high-statistics data samples, the ALICE experiment will identify, at central rapidities, J/y 
from b-decays, via the measurement of the pseudo-proper decay time distributions P31 . and will also 
reconstruct the Xc ^ J/W decay f46\. At forward rapidity, the contribution from b-decays will be 
estimated from the beauty cross section measurement carried out in the semi-leptonic decay channel [47.1 . 

7 Conclusions 

The ALICE experiment has measured inclusive J/ i/a production in the rapidity ranges \y\ < 0.9 and 2.5 < 
y < 4, through the decays J/y — > e+e^ and J/y M^M > respectively. The /Jx-integrated cross sections, 
based on data samples corresponding to integrated luminosities Li^ = 5.6 nb^' (for the J/y e+e^ 
channel) and Lim = 15.6 nb^^ (for J/y/ M^M") are aj/-^{\y\ < 0.9) = 12.4 ± 1.1 (stat.) ± 1.8 (syst.) 



Table 2: Summary of the results on the J/i/a differential cross sections. 



Pt Nj/ij/ Axe d^aj/yif/dp-rdy Systematic errors 



(GeV/c) 






(/ib/(GeV/c)) 


Correl. 

(^ib/(GeV/c)) 


Non-correl. 
(/ib/(GeV/c)) 


Polariz., CS 
(/ib/(GeV/c)) 


Polariz., HE 
(/ib/(GeV/c)) 


M<0.9 


[0;i] 


50±17 


0.122 


0.68±0.24 


0.02 


0.21 


+0.16,-0.18 


+0.08,-0.12 


[i;2] 


86±17 


0.076 


1.87±0.37 


0.07 


0.31 


+0.42,-0.50 


+0.28,-0.39 


[2; 3] 


79±13 


0.069 


1.89±0.31 


0.08 


0.23 


+0.33,-0.43 


+0.35,-0.44 


[3;5] 


75±13 


0.086 


0.72±0.13 


0.02 


0.09 


+0.06,-0.08 


+0.16,-0.13 


[5;7] 


50±9 


0.104 


0.40±0.07 


0.01 


0.05 


+0.001,-0.005 


+0.06,-0.08 


2.5 <3;<4 


[0;i] 


229±29 


0.280 


0.67±0.08 


0.06 


0.05 


+0.14,-0.21 


+0.13,-0.19 


[i;2] 


453±40 


0.287 


1.30±0.11 


0.12 


0.10 


+0.31,-0.35 


+0.19,-0.29 


[2; 3] 


324±26 


0.289 


0.92±0.07 


0.09 


0.07 


+0.17,-0.26 


+0.09,-0.19 


[3;4] 


253±21 


0.312 


0.67±0.06 


0.06 


0.05 


+0.12,-0.18 


+0.07,-0.11 


[4; 5] 


120±17 


0.359 


0.28±0.04 


0.03 


0.02 


+0.04,-0.05 


+0.02,-0.03 


[5;6] 


86±12 


0.392 


0.18±0.02 


0.02 


0.01 


+0.01,-0.03 


+0.01,-0.02 


[6; 8] 


80±12 


0.452 


0.07 ±0.01 


0.01 


0.01 


+0.007,-0.003 


+0.007,-0.008 


y 






daj/^/dj (jUb) 


m 


(Mb) 


(Mb) 


(/ib) 


[-0.9; 0.9] 


352±32 


0.085 


6.90±0.62 


0.28 


0.96 


+0.9,-1.3 


+1.0,-1.5 


[2.5;2.8] 


272±28 


0.117 


5.12±0.77 


0.49 


0.38 


+ 1.29,-1.62 


+0.94,-1.29 


[2.8;3.1] 


326±30 


0.383 


4.34±0.34 


0.41 


0.32 


+0.97,-1.06 


+0.85,-0.99 


[3.1;3.4] 


409±32 


0.469 


4.64±0.35 


0.44 


0.35 


+0.53,-0.92 


+0.46,-0.87 


[3.4; 3.7] 


271 ±26 


0.417 


3.59±0.30 


0.34 


0.27 


+0.57,-0.77 


+0.22,-0.53 


[3.7; 4.0] 


172±23 


0.215 


3.05±0.40 


0.29 


0.23 


+0.67,-1.01 


+0.09,-0.46 
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+ 1.8 (Ahe = 1) - 2.7 (Are = -1) Aib and c7j/^(2.5 < y < 4) = 6.31 ± 0.25 (stat.) ± 0.76 (syst.) 
+ 0.95 (Acs = 1) — 1-96 (Acs = —1) M^- The transverse momentum distribution was measured at 
both central and forward rapidity. Taking together the results from the muon and electron chaimels, the 
ALICE measurement of the inclusive J/y production cross section is particularly relevant in the context 
of charmonium studies at the LHC, for its coverage of both central and forward rapidities and for the 
lowest pt reach at y = 0. 
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